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We investigate theoretically how single- molecule spectroscopy techniques can be used to perform 
fast and high resolution displacement detection and manipulation of nanomechanical oscillators, 
such as singly clamped carbon nanotubes. We analyze the possibility of real time displacement 
detection by the luminescence signal and of displacement fluctuations by the degree of second order 
coherence. Estimates of the electro-mechanical coupling constant indicate that intriguing regimes 
of strong back-action between the two-level system of a molecule and the oscillator can be realized. 

PACS numbers: 07.10. Cm, 37.10. De, 42. 50. Hz 



Introduction Perfecting displacement detectors is the 
central motivation driving the research in nanomechanics 
[1] and a necessary step towards quantum manipulation 
of mechanical degrees of freedom. The main strategy 
is to couple the mechanical oscillator to a well controlled 
quantum system whose state can then be very accurately 
measured. A non exhaustive list of such detection sys- 
tems includes single-electron transistors [2], SQUIDs [3], 
Qbits [4], point contacts [5], optical [6], and microwave 
cavities [7]. A wide range of possibilities is opened by 
the emergence of a new class of detectors, where the dis- 
placement of the nano-oscillator modulates the energy 
splitting of a two- level system (TLS) [8, 9], which is then 
measured via optical resonance spectroscopy. In Ref. [8] 
the displacement (x) of the nano oscillator is detected by 
exploiting a Zeeman split TLS (Nitrogen Vacancy cen- 
ter) embedded on the tip of the oscillator moving in a 
strongly inhomogeneous magnetic field. 

In this Letter we propose using the single-molecule 
spectroscopy (SMS) [10-12] to detect the displacement 
of nano-oscillator s, such as carbon nanotubes suspended 
from the tip of an atomic force microscope, see Fig. 1. 
At liquid helium temperatures, the zero-phonon lines of 
single-molecule fluorescence excitation spectra are ex- 
tremely narrow, since the dephasing of optical transi- 
tion dipole due to phonons vanishes. For well chosen 
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FIG. 1: (Color online) Schematic representation of the pro- 
posed measurements setup: a) the geometry, b) modulation 
of molecular level splitting during the nanotube oscillations. 



fiuorophores/matrix systems [10], such as the dibenzo- 
anthanthrene DBATT embedded in a n-Hexadecane Sh- 
pol'skii matrix, the zero-phonon has a spectral width 
limited by the lifetime of the molecule excited state 
(- 10 - 20 MHz) [11, 12]. Under external electric field 
centrosymmetric molecules such as DBATT usually gain 
permanent dipole moments due to distortions induced 
by the surrounding solid matrix. This leads to a lin- 
ear contribution to the Stark shift, which is usually 
much stronger than the quadratic contribution. In disor- 
dered matrices such as in polymers the permanent elec- 
tric dipole moment can be as large as ID, and is around 
0.3D in n-Hexadecane Shpol'skii matrix (corresponding 
to ~ 3 MHz/(kV/m)) [12]. This allows one to use single 
molecules as highly sensitive probes of their nanoenvi- 
ronment and of local electric fields [13]. 

By setting a bias voltage between the nanotube (NT) 
and the conducting (and transparent) substrate, it is pos- 
sible to generate an electric field between the NT tip and 
the substrate at the limit of the discharge field Ec = 10^ 
V/m [14], with a very large radial gradient of the or- 
der of Ec/R = 10^^ V/m^, where ^ 1 nm is the NT 
radius [15]. Any small displacement of the NT will thus 
induce a large modulation of the molecular Stark shift, al- 
lowing an efficient mechanical/optical transduction. For 
typical molecules used in SMS experiments this gives an 
expected electric coupling constant a ~ 10^^ Hz/m, four 
orders of magnitude larger than the magnetic coupling 
constant observed in Ref. [8]. 

In order to show the capabilities of this detection tech- 
nique we provide explicit predictions for the luminescence 
spectrum and the second order photon correlation func- 
tion [g'''^\t)] of the molecule coupled to the oscillator. 
Specifically we show: (i) The modification of the molec- 
ular fluorescent excitation spectrum, (ii) The possibility 
of real time displacement detection, (iii) The relation be- 
tween g^'^\t) and the oscillator spectrum, allowing to use 
the detection scheme in the non-adiabatic regime, (iv) 
The possibility of the back action cooling of the oscilla- 
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tor. 

Model Our description of the proposed experiment is 
based on the Hamiltonian H = Hm + Hl + Hp^ot + 
Hose + H^. Here Hm = huj{x)n describes the electronic 
state of the molecule, modeled as a TLS, {|1), |2)}, with 
TT = |1)(2|, and h = |2)(2|. The TLS energy splitting is 
u{x) = ujQ -\- AjU ■ E{x)/h ~ cjo + A/i • E{0)/h + ax, 
where ujq is the bare level splitting of the molecule, 
A/i = (2|/i|2) — is its permanent electric dipole mo- 
ment, a = dx[Afl - E{x)]\x^o/^^ and E{x) is the electric 
field at the molecule that depends on the displacement 
X of the nanotube tip from its equilibrium position. For 
simplicity we will focus on a single oscillator mode, thus 
X is a scalar. The term Hl = l^L^re*^^^ + h.c describes 
the coupling of the TLS to a laser field of frequency ujl 
and with amplitude parametrized by the Rabi frequency 
The coupling with the photon vacuum environment 
is given by Hp^ot and leads to a decay rate F of the TLS. 
Finally Hqsc = +^^if^^/2 describe the 

selected mode of the mechanical oscillator (with resonat- 
ing frequency ujm and effective mass m) and its coupling 
to the a thermal environment of temperature T inducing 
a damping rate 7 and stochastic force fiuctuations. 

Under the usual assumptions of a short memory pho- 
ton vacuum [16] the equations of motion for the re- 
duced density matrix elements ai2 = (T21 = (tt) and 
(722 = 1 — cTii = (n) of the TLS take the form of Bloch 
equations: 



d(Tl2{t) 

dt 

dcr22{t) 
dt 



-i [6 + ax{t)] ai2{t) - -(Ti2{t) + 

i^^L [2(722 W-l], (la) 

-2QL^[ai2{t)]-ra22{t), (lb) 



with 6 = uj{0) — ul the detuning. 

On the same grounds, the dynamics of the average of 
X after tracing out the environment degrees of freedom 
is described by the Langevin equation: 



d^x{t) dx(t) 



dt^ 



+ 7- 



dt 



X{t) 



-^ha{t). (2) 



Here the force ^(t) is a gaussian fiuctuating field with 
{^{t)) = and {C{t)^{f)) = 2mjkBTSD{t-f) {ks is the 
Boltzmann constant and S]j{t) is the Dirac delta func- 
tion). We restrict to a classical description of the oscilla- 
tor {ksT ^ Tioom) that is the more relevant case for the 
typical experimental conditions. The last term in Eq. (2) 
describes the back-action force acting on the oscillator: 
mfba{t) = haa22{t)' 

We begin by assuming that the back-action force is 
negligible. The validity of this assumption will be dis- 
cussed later. 

Luminescence Spectrum Resonant fiuorescent excita- 
tion spectrum is proportional to the population of the 
TLS excited state (722- Solving Eqs. (1) to order ftj^ for 




FIG. 2: (Color online) Luminescence spectra for the TLS in- 
teracting with the oscillator for = 1-5 and ojm = 4r. 



r^i,/F <C 1 we obtain 



(3) 



Eq. (3) has to be averaged {{...) x) over the displace- 
ment fiuctuations describled by Eq. (2) to obtain the 
stationary population [15]. For small oscillator damping, 
7 <C cjM,r, this gives 



r + r. 
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(5 + nwM) +i(^ + ^<^)• 



:,(4) 



where Ik{z) are the modified Bessel functions, ("^ = 
o?{x?)xl^\/[ is an effective TLS-oscillator dimensionless 
coupling constant involving the temperature {{x'^)x = 
ksT /muo\^) and F^ = 27^^^ is the TLS mechanically- 
induced dephasing rate. 

Fig. 2 shows a typical luminescence spectrum obtained 
from Eq. 4. If not otherwise stated, in the Figures and 
estimates below we choose the following values of the 
parameters, F/(27r) = lOMHz, VLl = -ir, m = IQ-^^kg, 
and the mechanical oscillator quality factor Q = ujm h — 
10^. The main feature of the spectrum is the appearance 
of the sideband peaks at frequencies ^ ± nuju^ with n 
integer. Their intensity is controlled by through the 
Bessel functions. For <C 1 a measure of the oscillator 
fiuctuation can be obtained directly from the ratio of 
the height of two subsequent peaks /n+i (C^)/^n (C^) = 
C^/2(n + 1). Noteworthy that for low quality factors (cf. 
Q = 10 shown in Fig. 2) the peaks significantly broaden 
due to the mechanically induced dephasing [15]. 

It is instructive to compare this result to the case of 
a molecule driven by a perfectly periodic electric signal, 
such as a radio frequency (rf) field [11, 17]. Replacing 
OLx(t) with (rf^M cos{ujMt) in Eq. (3), where (rf is the 
dimensionless strength of the rf field, one obtains: 
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The overhne indicates averaging over the phase of the rf 
field. For ^.f <^ 1 the spectra due to Eqs. (5) and (4) 
essentiaUy coincide for = CV2. However, for (^rf ^ 1, 
as a function of n J^(Crf) oscihates, which results in 
the maximum of intensity shifting to the satellite peaks, 
whereas In{C^f/2) monotonically decreases. A simple in- 
terpretation of this fact can be obtained by averaging 
Eq. (5) over a gaussian distribution of [15]. One finds 
that the result is very similar to Eq. (4), apart from the 
broadening of the peaks. The spectrum (4) could thus 
be viewed with good accuracy as the ensemble average 
of many rf sources with random phases and a gaussian 
distribution of the squared amplitude. 

Real Time Position Measurement In the adiabatic 
limit, ujm <^ r, the sideband peaks merge into a single 
peak, whose maximum shifts proportionally to the oscil- 
lator displacement. Treating ax{t) as a time- independent 
shift in the Bloch equations, one obtains 



cr22{t) = 



ill 



[5 + ax{t)f + ^ 



+ 202 



(6) 



It is thus possible to monitor the real time position of 
the oscillator by measuring the variation of the lumines- 
cence intensity for a given value of the laser frequency: 
I[x{t)] = rjTa22{t), where 77 is the quantum efficiency of 
the detector. The change in the intensity in response to 
the oscillator displacement is 



G 



dl[x{t)] 



dx{t) 
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(7) 



Important characteristic of a displacement detector is 
its "sensitivity," i.e. the output noise referred back to the 
input parameter [1]: Sxx = \/ Sii{^) /G^ ^ where S'//(0) 
is the spectrum of the intensity ffuctuations at zero fre- 
quency. In the adiabatic regime one can treat these fluc- 
tuations as a Poissonian process of emitting single pho- 
tons, for which ^//(O) = I[x = 0]. We thus obtain 
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This expression takes minimum value at \5\ = VtL\/2 = 
r/2. Taking a = lO^^Hz/m and r] = 0.01 we ob- 
tain = ^/27/8^/T/r]/a - lOfm/VSz. This is 
on par with the best values of sensitivity obtained by 
non-interferometric displacement detection methods [18]. 
With respect to the latter, the important advantage of 
the proposed setup is its "open" nature, which allows 
a simpler integration of the mechanical resonator with 
others systems. 

Another essential characteristic of a detector is its res- 
olution. Ax = Sxx/"/^, 'i-^- the smallest displacement 
that can be measured [1]. Here At is the measurement 
time, which in the case of a damped harmonic oscilla- 
tor, Eq. (2), can be as long as 7"^ = Q/ujm- Thus, for 



the typical nanotube oscillator frequencies in the range of 
IkHz-lOMHz and quality factor Q = 10^ the resolution 
will vary from lOfm to 1pm, which again matches the 
values obtained by state of the art non-interferometric 
detection devices [1] . The device can be also used to per- 
form real time dispacement detection, without assuming 
a slowly damped oscillatory behavior of x{t). In order 
to have a sufficient signal to noise ratio in this case the 
limitation on the oscillator frequency is even more severe 
giving ujm ^ \^ / {olSxt)]^ = IMHz, where we used also 
the condition of validity of the linear expansion ax <^r. 

Second order photon correlation function In the most 
common situations, where one is interested in statis- 
tical characteristics, rather than the exact time evolu- 
tion of the oscillator, one can resort to the measurement 
of second order photon correlation function [^^^^(t,t')] 
[16], that is obtained by averaging over many two-photon 
events. We will show that this allows to improve the fre- 
quency resolution even beyond F. The measured signal 
for g^'^\t,t') reads: 



((7rt(f)7rt(t)7r(t)^(f)))^ 



(9) 



At equal times g^'^\t^ t) vanishes (tt^ = 0), since two pho- 
tons cannot be emitted by a single molecule simultane- 
ously (anti-bunching) . For large times {r = t — t' ^V~^) 
the photon emission events are no more quantum corre- 
lated and the degree of second order coherence reduces 
to the classical correlation function for the excited state 
occupation number: g'^'^\t,t') ^ {cr 22 (t) a 22 {t'))x/ {0-22)1- 
For vanishing coupling to the oscillator g^'^\t^t') and at 
lowest order in exhibits decaying oscillations at the 
detuning frequency 5: 



-Tt 



2 cos {6t) e 



(10) 



It is possible to obtain analytical expressions for 
g^'^\t^t') in presence of the oscillator at lowest order in 
^l/E < 1 and C < 1 [15]. The result is in general 
complex, but in some regimes it has simple and inter- 
esting interpretation. For 6 = ujm we find that the 
free oscillations described by Eq. (10) are in counter- 
phase with those induced by the oscillator leading to a 
reduction of the amplitude of the oscillations for tF ~ 1 
(cfr. Fig. 3). For Fr 1 the Rabi oscillations of Eq. 
(10) die out and the remaining time dependence is all 
due to the time-correlations induced by the oscillator 
{o22{t)o-22{t'))x- This is particularly interesting in the 
case of fast oscillators ojm ^ F 7 for which we find 



r2 



e ^2" cos (wmt - Q ^) . (11) 



The amplitude of g^"^^ gives thus a direct measure of the 
average amplitude fluctuation of the oscillator. In the 
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FIG. 3: (Color online) Top: g^'^\r) as a function of time 
for 6 = (jOm = 2r without /with coupling to the oscillator. 
Bottom: colorplot of g^'^\r) as a function of time and detun- 
ing. The dash line markes S = ujm- Parameters: ujm = 2r, 
C = 0.15, Q = 10^ 



opposite limit, ujm <^ T, one obtains 

9^'Ht,t') = l+(jy {x{t)x{t'))^, (12) 

valid for an arbitrary oscillator correlation function. 

Current- induced dissipation The strong coupling be- 
tween the NT and the TLS is possible due to the large 
electric field between the NT tip and the conducting sub- 
strate [15]. Oscillations induce thus dissipative electric 
currents. One can single out two different contributions 
from the currents flowing in the nanotube and in the 
substrate. We find that the dominant effect comes from 
the surface losses and that for typical parameters the in- 
crease in the oscillator damping rate is 1/ {z^ a sub) GHz 
nm^ S/m, where a sub is the substrate conductivity and z 
the distance from the substrate. Thus even if it is possi- 
ble to minimize this contribution by using higly conduct- 
ing substrates, in common situation it is likely that the 
quality factor of the oscillator will be reduced. 

Back-action and the effective temperature Let's now 
discuss in more details the role of the back-action force. 
Its effect is controlled by two parameters, the intensity 
of the force (measured by the ratio between the displace- 
ment induced in the equilibrium position of the oscillator 
ha/mujj^ and \/ {x'^)x), and the frequency and duration 
of its fluctuation. In order to estimate its average effect 
taking both parameters into account we resort to the lin- 
ear back-action theory [19]. Within this approach action 
of the TLS on the oscillator can be obtained by evalu- 
ating the fluctuation spectrum of the force {hah) in the 
absence of the mechanical oscillator. This is proportional 
to 



SnnH = j dte'^' (n(t)n(O)) 



r(n) 



((5-0;)" + 



(13) 



where (n) = ^1/ {5^ ^T^ /A), h 



(n). For Q > 1 



be fully characterized by an effective temperature (T^. 
deflned by 



coth 



/ hUM 



7 coth (^)+ 71 coth (,^^) 



7 + 7i 



one can show [19] that the oscillator fluctuations can then 



^(14) 

where Ttls = hunS:^^/ [2kBS-^], 71 = ( ax^pf)^ S'^ , 

^tn Snni^M) ^ Snn{-^M), and X^pf ^J h / {2mUJ m) - 

When 7i ^ 7 the oscillator temperature is completely 
determined by TtlSi which in turn is controlled by the 
detuning, 5. For positive values of (5, i.e. for the pho- 
ton energies below the transition energy, this results 
in cooling the oscillator down to the temperature Ttls 
[20, 21]. At the optimal value 5 = y^uj]^ + /4 we 
obtain Ttls = ^F/(4/cb) for ujm <C F and Ttls ^ 
^^m/[(2/cb log(4co'M/F)] in the resolved side-band limit, 
ujm ^ F. Note that unlike in the case of coupling to 
a photon cavity [20, 21], where one can inflnitely in- 
crease the number of cavity photons, the population of 
the excited state ((n)) entering linearly the expression for 
7i saturates at high laser strengths and remains smaller 
than one half. 

The effective temperature expression (14) holds for 
7i <C F [20], whereas cooling requires that 71 7. In 
the resolved side-band limit these conditions result in [15] 



{n)^ whereas for ujm <^ T 
the conditions read the same, but with x^pf replaced by 
y^h/{2mr). For ujm ~ F we thus obtain that the lin- 
ear response description is correct for a < lO^^Hz/m, 
i.e. under all realistic circumstances, whereas the back- 
action effects are negligible for a <C lO^^Hz/m. This 
gives the quantitative condition for discarding /ba in the 
flrst part of this paper. This also indicates that there 
is a window of at least two orders of magnitude of the 
coupling constant where back-action can be used to ma- 
nipulate conveniently the oscillator. The dependence of 
the effective temperature on 6/r and cjm/F is shown in 
Fig. 4: The temperature of the oscillator can be reduced 
by a factor of a 1000 when cjm/F ~ 0.1 — 0.4. 

Conclusion In this Letter we have shown that single- 
molecule spectroscopy can be used as an efficient tool for 
detecting and manipulating nano-oscillators, such as car- 
bon nanotubes attached to an AFM tip and interacting 
with the molecule via electrostatic forces. We analyzed 
different detection schemes which allow obtaining either 
real time or averaged information on the oscillator dis- 
placement. Estimated very large values of the coupling 
constant allow reaching high detection efficiency. Fur- 
thermore, the back-action of the molecule on the oscilla- 
tor can be made sufficiently large to enable strong cooling 
of the oscillator with the proposed set-up. We considered 
in this Letter only the linear part of the TLS-oscillator 
coupling, but by properly position the nanotube it is pos- 
sible to work in a regime where the dominant coupling 
is quadratic, leading to the wide range of possibilities in- 
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FIG. 4: (Color online) Ratio kBTosc/ {fiojM) as the function of 
the detuning and the oscillator frequency for a — 10^^ Hz /m, 
T=1K ^ 22GHz. 

eluding that of measuring the energy of the oscillator [22] . 
Single-molecule spectroscopy appears thus as a promising 
new technique for detecting displacement and exploring 
strong coupling regimes of a driven two-level system with 
a nano- mechanical oscillator. 
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We provide here more detailed derivations of some of the results presented in the main manuscript^ 
PACS numbers: 



I. SPECTRUM OF THE TLS COUPLED TO AN OSCILLATOR 



Solving Eqs. (1) and (2) of Ref. 1 to order Ql for Ql/T ^ 1 (and neglecting the back-action) we obtain 

J —OO 



(722 (i) = -20iS 



r rfiiai2(ti)e-^(*-*^) 

J —OO 



(la) 



(lb) 



Substituting Eq. (la) into Eq. (lb) and using algebraic manipulations one can show that (722 (^) is given by Eq. (3) 
of Ref. 1, which is convenient for the analysis of the adiabatic limit and the degree of second order coherence. For 
now we focus on Eqs. (1). One can use the Gaussian properties of the stochastic noise, ^(t), (and therefore also of 
x{t)) to carry out the averaging over the oscillator position fluctuations: 

^^-^^ St\ drx{r) ^^ ^ e^''^~^^\ (2) 

(3) 
(4) 



where we took into account that {x{t)) = 0, whereas 

$(*) = -— dT2 dn{x{T2)x{n)), 
with the oscillator correlation function given by 



{x{t)x{t')) 



{x'}e- 



-f\t-t'\ ( 

' I 



cos[co'm(^ — t')] + 773 — sin[u;M|^ — 
2ujm 



Here A = —7/2 + ic^M, = a/'^m ~ 7^/4, and {x^) = ksT / {muj\^) is the level of the oscillator fluctuations as 
determined by the equipartition theorem. 
We now have 



{(^22) 



2nie-(-+-*) ^ 



dte 



i[(5-j(r+r^)]t+;^e^*+;^*e^** 



(5) 



where = 27^^ is the dephasing rate of the TLS induced by the mechanical noise, z = C^(A*)^/(2ia;M^M), and 
= a^{x^) /uj\j- is the dimensionless coupling strength between the oscillator and the TLS. 
Expanding the integrand of Eq. (5) in series in ze^^ ^ z'^e^ * one can carry out the integration with the result 



21^?e-(^+^*) 



2" {z*y 



n\m\ i[S^{m- n)cuM] + ^ [F + F^ + (n + 771)7] 



This equation was used to obtain Fig. 2 of Ref. 1. 

In the limit of the high oscillator quality factor, Q = ujm/i ^ 1, we have z ^ and A 



cos(cc;m^), and can use the expansion in terms of the modified Bessel functions^ 



COS 



ike 



(6) 

CUM • Then 
(7) 



k= — oo 



to obtain Eq. (4) of Ref. 1. Comparing with Eq. (6) we see that this corresponds to assuming that 7 <C com and 
(n + m)7 <C F + F^, i.e. the approximation is valid when one can neglect the satellites with indices higher than 

(r + r^)/7- 
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II. CORRESPONDENCE BETWEEN OSCILLATOR AND AC SPECTRA 

Let us consider a molecule driven by an undamped oscillator with a constant amplitude A: x{t) = A cos [uJMt + 
The resulting spectrum is given by Eq. (6) of Ref.l, where the argument of the Bessel function is (rf = Q^A/cjm- We 
now will average this spectrum over the thermal distribution of the oscillator energies 



We{E) 



1 



(8) 



The oscillator energy is related to its amplitude as E{A) = muj\jA^ which allows us to write the probability 
distribution for the amplitude 



Wa{A) 



dE{A) 
dA 



We[E{A)] 



(Note that (A^)^ = 2(x2)^.) 
We now evaluate the average 



aA 



f 

Jo 



dAWA {A) 



(9) 



(10) 



where {Crf)x = <^^(^^)x/^m5 and we have used the relation'' 



1 



af3 



/ dxe ^^'^^xJp{ax)Jp{/3x) = -^e \ o 2 

Jo 

The resulting averaged spectrum coincides with the one given by Eq. (5) of Ref. 1 (at zero damping = 0). 



(11) 



III. DEGREE OF SECOND ORDER COHERENCE 

A. Calculating the degree of second order coherence, g'^'^\r) 

The usual way to obtain g^'^^r) [Eq. (10) of Ref. 1)] is by using the quantum regression theorem, ^'^ which relates 
g^'^\r) to the evolution of the average dipole moment and the excited state occupation number. This is done by 
deriving equations for the correlation functions 

CMt,t') = {7r\t')n{t)Tr{t')), C22{t,t') = {7r\t')n\t)n{t)Tr{t')), (12) 

using the same approximations as those that were used in deriving the Bloch equations [Eqs. (1) in Ref. 1]: 

(13a) 



lCi2(t,i') = -i[S + ax{t)]Ci2{t,t') - ^Ci2{t,t')+inL [2C22{t,t')-a22{t')]. 



d_ 
Ft 



C22M = -2f]z,S[Ci2(t,i')]-rC22(i,t')- 



(13b) 



Eqs. (13) are solved in the interval \t' ^t] with the initial conditions Ci2{t' ^t') = 0, and C22{t'^t') = which readily 
follow from the definition of these functions. 

So far we have not performed the statistical averaging over the thermal fluctuations of the oscillator, x{t). Denoting 
such averaging by {...)x we define the degree of second order coherence as 



{C22{t,t'))^ 
{'T22{t))^{<J22{t'))^' 



(14) 



Implied in this definition is that the Brown- Twiss experiment does not directly measure g^'^\t, t'), but rather measures 
separately C22{t, t') and (722{t)^ which are averaged over the oscillator fluctuations after many repeated measurements, 
and later operated numerically to form g^'^\t^t'). Note that after the averaging the degree of second order coherence 
will depend only on the difference of the two times: t = t — t' . 
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B. Degree of second order coherence to order 



When solving Eqs. (13) to order C22{t,t') is given by 



(15) 



where (J22{t) is given by Eq. (3) of Ref. 1. Thus, C22(^,^0 actually of order Vt\^ but the lowest order contribution 
to the degree of second order coherence, defined by Eq. (14), will be of zero-th order in ^l^- 

For a Gaussian statistics of x{t) one can explicitly carry out the statistical averaging of (J22{f) and C22(^, t')- In the 
following we limit ourselves to the limit C <C 1 to order where most of the integrals can be evaluated analytically. 
The expressions given below were used to obtain Fig. 3 of Ref. 1. 

After straightforward but tedious calculations one arrives at 



{(J22{t)(J22{t'))^ 
{(^22)1 



{(^22Y 



(^22)^ 



(^22)^ 



1 + -23? 
1 + ^23? 

+ ^ 



2 



2/q 



1 + 23? 



~2 ^0 



^-f r 



(16a) 

(16b) 
(16c) 



(<T22(i)<T22(i')). + (<Ti2>.e-r(*-*') 



where we have defined 



(^22)^ 



cos {b(t - 



(^22(^)^22(^0), 
(^22)^ 



5R 



i (/+ + /_ - 2/0) 



2(^ 



sin [^(t - t')\ ! 



" 2 



+ 



K/+ + /--2/0) , 2((5+f) 



Jo 

r+00 

Jo 

nt-t' 

= dT{x{T)x{0))^, 
Jo 

-t') = -y dT2 J^^ dri {x{T2)x{n))^ . 



(17) 

(18a) 
(18b) 
(18c) 
(18d) 



Note that for t — = we have I± = /q, and $(0) = 0, so that the coherence function takes zero value, which 
corresponds to antibunching. 



C. Vanishing coupling to the oscillator 

When the coupling to the oscillator vanishes the degree of second order coherence takes the simple form 

^(2) (r) = 1 + e-^^ - 2 cos (Sr) e' (19) 

which illustrates all the main features of ^^^Hr) originating from single molecule luminescence: i) antibunching, i.e. 
^(2)(o) = 0; ii) unitary limit at long times, g^'^\r -\-oo) 1; and in) oscillations with the frequency of detuning at 
times smaller than F"-^. These are illustrated in Fig. 1 and Fig. 2 (left panel). 
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FIG. 1: Degree of second order coherence, g^'^^r), for a TLS not coupled to an oscillator as a function of time for various values 
of detuning. 



D. Adiabatic limit 



Let us now consider the adiabatic limit, i.e. when the correlation function {x(t)x(t'))x changes on a time scale much 
slower than . In this case the correlation function can be taken out of the integrals of Eqs. (18) and, for times 
larger than the lifetime of the excited state of the molecule, r ^ F"^, we obtain Eq. (11) of Ref. 1. 

In the adiabatic limit it is possible to obtain more general results, not restricted to small values of the effective 
coupling Firstly, for the time intervals t — t' > F"^ relevant to the adiabatic limit, the degree of second order 
coherence is reduced to the correlation function for the excited state population of the molecule: 

where a22 is given by Eq. (7) of Ref. 1. One can carry out the statistical averaging and express the result in terms 
of a double integral: 

(^22(^)^22(0).= f^j J dr2j dricos[(5(r2+ri)]e-^-^ ^ M^MO)r2r.^ (21a) 

q2 r-\-oo ^ 2/ 2x 2 

(^22^). = -ff / dTCOs{ST)e-^-^^^. (21b) 

^ J — OO 

Eqs. (21) show that in the adiabatic limit g^'^^r) is a non-linear transformation of the oscillator correlation function, 
{x{t)x{t')). Thus, although g^'^\r) may not always correctly reproduce the oscillator correlation function, it allows to 
deduce characteristic frequencies of the oscillator. 

In the limit where the TLS-oscillator coupling is stronger than the molecular linewidth, a^{x'^) ^ F^, Eqs. (21) 
can be simplified by neglecting the terms F|r|/2 in the exponents. In this case the intergrals can be evaluated exactly 
and we arrive at 



g^^\t,t') = =exp 

0)' 



/ {xit)xit') 



{x{t)x{t')) 



a^{x^) {x^) ^ {x{t)x{t')) 



(22) 



whereas 



which is the population of the excited state of the molecule with inhomogeneously broadened line. 
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E. Fast oscillator 



At times t — t'^T ^ the degree of second order coherence, Eq. (17), is reduced to the correlation function for the 
population of the excited state of the TLS: 



Ml 



1 



+ 00 



JZ 72 COS [Uj{t -t)], 



27r 



where we have re-expressed the integrals of Eqs. (18) in terms of the oscillator spectrum, defined as 



+ 00 



dte'-' (x(t)x(O)), . 



In case of a linear oscillator with correlation function given by Eq. (4) this spectrum is given by 



(w2-a;2^)^+wV' 



and we obtain 

(^22)1 



£2 

4 



(24) 



(25) 



(26) 



(27) 



In the limit 5 = ujm ^ T the second term in this equation can be omitted and we arrive at Eq. (13) of Ref. 1. 
For negligibly small oscillator damping Eq. (27) becomes 



(^22(^)^22(^0), 
Ml 



1 + C^duj 



M 



1 



r2 



r2 



cos [i^M{t- t')] = 1 + Acos [a;M(^ - ^0] • (28) 



If we now consider the degree of second order coherence at times where e ^'^/^ ^ we can omit the terms proportional 



to C^e ^'^/^ and e but keep the oscillatory term 2e ^ 2 ^ cos [5(t — t')]: 

gi2)(^t,t') = ("22(^)"2f^))- _ 2e-^ cos m - t')] . 



(29) 



When 5 = ujm the two terms oscillate at the same frequency but have different signs, i.e. the overall oscillation is 
suppressed: 
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(2) 



(r) = 1 + (^A - 2e"^) cos {5t) . 



The amplitude of the oscillations thus vanishes at 



To = - log 



log 



2q^2(^2\ 



(30) 



(31) 



Since zeros of the cosine function {ujm^ = 7r(n + 1/2)) occur at equally spaced intervals, tq is easily identified as the 
additional crossing of ^^^^ = 1. The exception is the case when tq coincides with one of the zeros of the cosine. Then 
the two zeros merge and the oscillations only touch ^^^^ = 1. (^^^^ = 1 can be identified as the asymptote of the 
oscillations at times r ^ F"-*^.) 



IV. ESTIMATE OF THE ELECTRIC FIELD GRADIENT 



In Ref. 1 we estimated that the electric field gradient due to the nanotube by taking the ratio of the maximum 
electric field (discharge field, Ec = lO^V/m) to the nanotube radius {R ~ Inm): Ec/R = lO^^V/m^. Here we provide 
more rigorous electrostatic calculations which support this simple estimate. 

We treat the nanotube as a perfectly conducting cylinder of radius R with infinitely thin walls, suspended above a 
conducting surface (substrate). The boundary of the substrate is specified by 2; = 0. The nanotube axis is taken to 
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FIG. 2: Degree of second order coherence, g^'^\r)^ as a function of time and the detuning for a molecule not coupled to an 
oscillator (left) and a molecule coupled to an oscillator with = .1, um — 2r, and Q — 10^ (right). 



coincide with the z-axis, its length is L, and the distance between the bottom end of the nanotube and the surface is 
zq. The nanotube has potential V in respect to the surface, while the surface itself is assumed to have zero potential. 
The electric potential due to the charge in the nanotube and the image charges induced in the substrate is^ 



(^(r,z)= dz'\{z')[F{r,z-z')-F{r,z^z')], 



where X{z) is the linear charge density along the nanotube (we use cgs units), and the integration kernel is 



F{r,z)= / 
Jo 



dcj) 



1 



27r j^R2 _ 2rR cos 0-t Z^ 



(32) 



(33) 



Since we assume that the nanotube is perfectly conducting, all its points have the same potential, i.e. 

(p{R, z) = V for zo < z < zo^ L, (34) 

which is the integral equation for determining the distribution of the electric charge on the nanotube. Note that Eqs. 
(32) and (34) scale linearly with the potential difference, i.e. increasing/reducing V will proportionally change the 
electric field affecting the molecule. 

The magnitude of the electric field at the surface, obtained by numerically solving Eq. (34), is shown in Fig. 3 (left 
panel) as a function of the radial distance, r, from the nanotube axis, and for different different values of the distance 
between the tip of the nanotube and the surface, zq. We used the following parameters: V = lOmV, R = Inm, and 
L = bOR. The value of the nanotube length that we used was found to be bigger than the intra-nanotube screening 
length, i.e. the results are valid also for much longer nanotubes. 

The corresponding radial gradient of the electric field is shown in the right panel of Fig. 3 and agrees by the 
order of magnitude with the simple estimate provided in Ref. 1. The coupling constant is obtain by multiplying this 
gradient by a typical molecular Stark shift (e.g. IMHz per every kV/m for dibenzanthanthrene (DBATT) molecules 
in hexadecane (HD)^.) 



V. DISSIPATION 



We now estimate the dissipation due to the Joule's heat originating from the currents induced in the the substrate 
and the nanotube when the latter vibrates. For simplicity we assume that all of the nanotube charge is accumulated 
in its tip, in the region with length of the order of the nanotube radius. This charge then can be estimated as^: 
Q{t) ~ X{0,t)R ~ VR/\og [4:z{t)/R]^ where z{t) is the distance between the tip and the substrate. 

We model the displacement of the nanotube tip as that of a straight rod in xz plane: x{t) = Xmax ^os (Qt), 
z{t) = zq -\- L — L? — x^{i) + x'^{t)/{2L). (We cannot use here the Euler-Bernoulli beam model, "'^^ since it 

ignores the displacement along the nanotube, i.e. effectively elongates it.) 

Dissipation by the intra-nanotube currents: The instantaneous power for the Joule heat generated by the currents 
flowing to and out of the nanotube tip is W{t) = pNTL[Q{t)\^ ^ whereas the average energy dissipated during the 
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period of the oscillations is 



Wtube 



Pnt 
8L 



.olog(^; 



2 4 



(35) 



where Qo = VR/\og [Azq/R], pnt is the nanotube resistance per unit length. 

Note that the dissipated power is proportional to the forth power of the oscillations amplitude, which corresponds 
to a non-linear dissipative force Fjoui = —^"iJouix^x (rather than the usual viscous friction force F^isc = —mjoc). 
This nonlinearity results from the geometric symmetry: the flow of charge in the nanotube is due to the change of the 
distance between the nanotube and the substrate, z{t) — zq, which is proportional to second power of the oscillator 
displacement, x{t). A linear term would appear, if the equilibrium position of the nanotube was not perfectly vertical. 

Dissipation in the substrate: The nanotube oscillations induce redistribution of screening charge in the substrate, 
i.e. a flow of currents, which also dissipate energy via Joule heat. For the typical nanotube oscillator frequencies 
(smaller than IGHz) the motion of the nanotube and screening charges is still sufficiently slow to be described as a 
quasi-electrostatic problem. The Poisson equation now contains time-dependent charge density and potential, and 
needs to be supplemented by the appropriate time-dependent boundary conditions (relating for the instantaneous 
electric field outside and the inside the substrate to the surface charge), the current continuity equation, and the 
consistency relation between the current density and the electric field: j = (Jsub^ {o'sub is the substrate conductivity.) 

The resulting solution can be shown to be similar to the one obtained by the method of images, but with the 
retarded image potential. The dissipated energy is then obtained by integrating the local density of the Joule power, 
w = jE = (JcE'^^ over the volume of the substrate. The energy dissipated in the substrate per a period of the nanotube 
oscillations is then 



W 



2^2 



substrate 



'ibQo'^ 



3zl 



M 



(36) 



10 13. 



The amplitude of thermal nanotube oscillations can be determined from the equipartition theorem: 
{x^) = ksTx /{mujj^). The frequency of the main nanotube mode can be taken from the Euler-Bernoulli theory 

/m = ^M/(27r) = ji^l^J^~{d^~^-~d^/{S7^L'^)^ where /ii = 1.875, E is the nanotube Young modulus, p = 2300kg/m^ is 

the density of graphite, d^, di are the outer and the inner diameters of the nanotube. We then obtain the estimates of 
the dissipation due to the intra-nanotube and substrate currents shown in Fig. 4, where we took R = Inm, E = 0.4 
TPa^^, Pnt = 6M0m/m, and Tk = IK. The quantity shown in Fig. 4 is the ratio of the oscillator energy to the 
energy lost per one period of the oscillations due to the given dissipation mechanism, i.e. the effective quality factor. 
We observe that the dissipation due to the substrate current, shown for several values of the substrate conductivity, 
is the dominant of the two dissipation mechanisms. It can be quantified in a more practical way by defining the 
damping rate associated with energy loss of Eq. (36): 



Isub 



127rmzncr 



Q^^sub 



(37) 



This equation is valid in the gaussian as well as in the SI units, and we took account that the conductivity of IS/m 
corresponds to frequency of 9GHz, much greater than the possible frequency of the nanotube oscillator. For the 
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FIG. 4: Fraction of the oscillator energy lost during one oscillation period due to the intra substrate currents (for ac — 
1, 10^, 1^^ S/m) and due to the currents generated within the nanotube. 

typical values of the parameters this equation gives the estimate ^suh = IGHz nm^ S/m, which means, e.g., that 
for /m = IMHz oscillator removed by = lOnm from the substrate surface we need the conductivity of at least 
(^suh — 10^ S/m to ensure the quality factor of Q = 10^. 

VI. WEAK COUPLING, STRONG COOLING, AND WEAK BACK-ACTION 

To order Vt\^ the fluctuation spectrum for the population of the excited TLS state is given by 



^nnM^ ^ die'"'' {n(i)nm 



^^-^ (38) 



((5-c^)" + 



where n — n — (n). Although it is possible to obtain the expression for this spectrum applicable to all orders in (7/,, it 
is rather cumbersome and would complicate comparison with the case when the oscillator is cooled by coupling to a 
photon cavity. ^^'^^ However, one should keep in mind that the average population of the excited state of the molecule 
saturates at high laser strengths: 



The effective oscillator temperature is given by the following expression 



15 



/ \ 7coth i^p^ + 71 coth f ^.^^^ 



where the cooling rate is determined by 

71 = [Snn{^M) - Snn{-^M)] = /^%f{n)V 

whereas the effective TLS temperature is defined by 

71 coth (^ 2kBTTLs ) ^ t^^^*^^^) + Snn{-^M)] = /^%f{n)T 



1 



1 1 



r2 



(41) 



(42) 



In the above expressions we have defined the coupling parameter A^^j = axzpf^ where Xzpf = ^/fl/ {2mujM) is the 
zero-point fluctuation length for a quantum oscillator of frequency ujm- 

We now specify three conditions: the (i) weak coupling condition, 71 <^ F, which specifies the limit where the 
effective temperature description by Eq. (40) is applicable^^; (ii) the condition of absence of cooling, 71 <C 7; and 
(iii) the condition that one can neglect the effect of TLS fluctuations on the oscillator, JiTtls <^ (which is the 
high-temperature limit of 71 coth [Tiujm / {^ksTTLs)] ^ 7 coth [Tiujm / {2kBT)].) 
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In the resolved side-band limit, uom ^ T, the first term in each of Eqs. (41), (42) takes maximum value at (5 = ujm^ 
whereas the other term is negligible and one arrives to the following conditions 

(43a) 
(43b) 

(iii) A,^^«./;^^^ = Jf^W^. (43c) 

Note that, although the excited state population (n) is also a function of the detuning, we treat it as an independent 
variable, since it can be adjusted by changing the laser strength. Yet, we stress that (n) cannot exceed one half, which 
means that (unlike in the case of cooling with a photon cavity) one cannot at will change the validity of conditions 
Eq. (43c). 

For the parameters used throughout our manuscript and for ujm = T = 10^ Hz we have x^pf ~ 20 pm (which is 
about two orders of magnitude smaller than ^/{x^ ~ Inm.) With the best possible value of = lO^^Hz/m this 
gives us Azpf ~ 10^ Hz. Given that (n) ^ 1 this means that we are always safely in the regime of weak coupling. 
H owever, as we see from the second condition, weak coupling does not exclude strong cooling^ i.e. the regime in which 
/ {n) ^ ^zpf ^ r/ a/ (n). Thus, the proposed setup can be used for cooling. 

In absence of cooling, i.e. when axzpf <C (^), the third condition is automatically satisfied, since 

V^^^ /^zpf ~ 50, which means that one can neglect the TLS back-action on the oscillator for values of about 
a = lO^^Hz/m and less. 

In the limit where the side-bands are not resolved, i.e. cum T, maximizing 71 and ^iTtls in respect to 5 results 
in conditions 

(i) Ar « (44a) 



(ii) Ar«W7\, (44b) 

V (") 

which differ from Eqs. (43c) by the replacement Xzpf xr = \/h/{2mT)^ ^zpf ^ ^r- Since in the previous estimate 
we took (jOm = r, we come to the same conclusions: one can neglect by the TLS back-action up to a = 10^''Hz/m, 
whereas for higher values one realizes regime of strong cooling, which however can be studied within the weak coupling 
approach employed here. 
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